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Abstract: Recent results from the WMAP measurements of the cosmic background radi- 
ation yield very tight constraints on the relic density of supersymmetric cold dark matter. 
We combine the WMAP constraint with those from the anomalous magnetic moment of the 
muon and the b — > 57 branching fraction in a x 2 determination over the minimal supergrav- 
ity model (mSUGRA) parameter space. The most favored region of mSUGRA parameter 
space for almost all tan/? values is the hyperbolic branch/focus point (HB/FP) region, 
with moderate to small values of superpotential Higgs mass |/x| and large GUT scale scalar 
mass rriQ. These favored regions of mSUGRA parameter space can be probed by direct 
search experiments for supersymmetric dark matter. An exception to the HB/FP region 
can occur at very large tan/3 with positive fi values, where wide regions allow resonance 
annihilation of neutralinos in the early universe. 
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The past decade has witnessed increasingly precise measurements of the anisotropics 
of the cosmic microwave background radiation left over from the Big Bangp]]. The most 
recent results come from the Wilkinson Microwave Anisotropy Probe (WMAP) satellite 
measurements. Astonishingly, an analysis of their results pinpoints the age of the universe 
to be 13.7 ± 0.2 GyrsQ. In addition, the geometry of the universe is flat, consistent with 
simple inflationary models. The dark energy content of the universe is found to be about 
73%, while the matter content is about 27%. A best fit of WMAP and other data sets to 
cosmological parameters in the AC DM cosmological model yields a determination of bary- 
onic matter density Vt b h 2 = 0.0224 ± 0.0009, a total matter density of VL m h 2 = 0.135^;^, 
and a very low density of hot dark matter (relic neutrinos). From these values the cold 
dark matter (CDM) density of FlcDAih 2 = 0AU6±%%\1\ (at 2a) can be inferred. The 
new WMAP results can thus be used to obtain more severe constraints on particle physics 
models that include candidates for cold dark matter, such as supersymmetric theories |J. 

It is well known that the lightest supersymmetric particle (LSP) of many supersym- 
metric models has the necessary attributes to make up the bulk of cold dark matter in 
the universe This holds true especially in supergravity models where supersymmetry 
breaking occurs in a hidden sector of the model||]. SUSY breaking is communicated to the 
observable sector via gravitational interactions, leading to soft SUSY breaking mass terms 
which can be of order ~ 1 TeV, so that the hierarchy between the weak scale and any other 
high scale such as Mqut or Mp\ can be stabilized. The simplest of these models assumes 
a flat Kahler metric and a simple form for the gauge kinetic function at the high scale. 
Here, motivated by gauge coupling unification at Mqut — 2 x 10 16 GeV, we thus assume 
common scalar masses mo, common gaugino masses mi/2, and common trilinear terms 
Aq all valid at scale Q = Mqut- Below Mqut-, the effective theory is assumed to be the 
minimal supersymmetric standard model (MSSM). The weak scale sparticle masses and 
couplings are determined by renormalization group running between Mqut and M wea k, 
which leads to radiative electroweak symmetry breaking (REWSB). The mSUGRA model 
parameters 

m o, m i/2) Ab tan/3 and sign(fi) (1) 

then determine all superparticle and Higgs boson masses and mixings. Here, tan/3 is as 
usual the ratio of Higgs field vevs. We use the program ISAJET V7.64P 1 for our sparticle 
mass calculations |J. 

Once the sparticle masses and mixings are determined, a variety of observable quan- 
tities can be calculated. In this letter, we focus especially on the neutralino relic density 
Qg^h 2 . The relic abundance of neutralinos can be calculated by solving the Boltzmann 
equation as formulated for a Friedmann-Robertson- Walker universe. Central to this cal- 
culation is the evaluation of the neutralino annihilation and co-annihilation cross sections, 
which must then be convoluted with the thermal distribution of neutralinos (and possibly 
other co-annihilating particles) present in the early universe. We adopt the calculation of 
Ref. Q, wherein all relevant annihilation and co-annihilation reactions are included along 
with relativistic thermal averaging || (see also Ref. || for a recent relic density calcu- 

1 Isajet 7.64p is Isajet 7.64 modified to gain access to low fj, sparticle mass solutions. 



- 1 - 



lation). Four regions of mSUGRA model parameter space emerge where the CDM relic 
density is consistent with measurements. These include A.) a bulk region at low mo and 
low mi/2; where neutralino annihilation occurs mainly via t-channel slepton exchange, B.) 
the stau co-annihilation region where Z\ —f\ and f\ — f\ annihilations contribute [10], C.) a 



region where 2m^ ~ mA,H, where neutralinos can annihilation via s-channel pseudoscalar 



(A) and heavy scalar (H) Higgs bosons[ll], and D.) the region at large itlq where 



becomes small [ 12 1 (known as the hyperbolic branch/focus point (HB/FP) region []l3|, H]), 
and the growing higgsino component of Z\ allows for efficient neutralino annihilation and 
co-annihilation [ 15 ] . 2 



Another important constraint on the mSUGRA model comes from comparison of the 
predicted rate for b — > 57 decay against experimental measurements. Here, we adopt 
the branching fraction BFib — ► S7) = (3.25 ± 0.54) x 10~ 4 for our analysis, and use the 
theoretical evaluation given in Ref. fl7|, Generally, the value of BF(b — * 57) calculated 
in the mSUGRA model differs most from the SM prediction in the region of low mo and 
J71W2, where sparticle masses are light, and SUSY loop contributions can be large. 

The recently improved measurement of the muon anomalous magnetic moment = 
{g — 2) M also provides an important constraint on supersymmetric models A recent 
determination of the deviation between the measured value of and the SM prediction 
has been made by Narison, including additional scalar meson loops |^]. His determination 
using e + e~ — > hadrons data to evaluate hadronic vacuum polarization contributions yields 
Ao^ = (24.1 ± 14.0) x 10 -10 , which we adopt for this analysis[21]. An alternative determi- 



nation using r-decay data may include additional systematic uncertainties, and is usually 
considered less reliable. 

Finally, we include in our determination of allowed parameter space direct superparticle 
and Higgs boson search results from the LEP2 experiments. The most important of these is 
that rriyy > 103.5 GeV on the lightest chargino[^2|, and mh > 114.1 GeV when the lightest 



Higgs boson is SM-like[23]. In regions where tua is small, this bound may be considerably 



reduced to rrih ~ 90 GeV, depending as well on the value of tan (3. 

In this analysis, we compute a x 2 value constructed from the mSUGRA model calcu- 
lated values of fJ^/i 2 , BF(b — > 57) and a^ USY , along with the above mentioned central 
values and error bars (24)]. If the relic density ^^/i 2 falls below the WMAP central value, 
then we do not include it in our \ 2 determination since other forms of CDM may be present. 
Thus, at each point in model parameter space, we determine the value of x 2 > and represent 
the value log(x 2 /3) by various colors in the itlq vs. m^j plane for different values of tan/3 
and sign(fi). The green regions generally correspond to a x 2 /dof value less than 4/3, while 
yellow regions have 4/3 ~ x 2 /dof ~ 25/3. The yellow regions shade into red for larger 
X 2 /dof values. We adopt = throughout our analysis. In general, our conclusions do 
not change qualitatively upon variation of Aq, unless extreme values of the parameter are 
adopted. 



2 In this paper, we consider mSUGRA solutions with Ao = only. For particular Aq values, the value 
of can be dialed to near degeneracy with , so that a fifth region of stop- neutralino co-annihilation 
occurs llfitl. 
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In Fig. [T], we show the mo vs. m-y/2 plane for < and values of tan/? = a.) 10, b.) 
30, c.) 45 and d.) 52. The gray regions are excluded either by a non-neutralino LSP, or 
by a breakdown in the REWSB mechanism, signaled by /u 2 < or by m\ < 0. The blue 
shaded region denotes points excluded by the LEP2 bound > 103.5 GeV, and the 

region below the blue contour is where m/j < 114.1 GeV. We see from frame a.) that only 
tiny green regions occur along the excluded region at low mo where stau co-annihilation 
occurs, or at the boundary of "No REWSB", where \fj,\ becomes small, the HB/FP region. 
The bulk region at low mo and mm is excluded by the LEP2 bound on m^, but also has a 
large negative value of a^ USY and a large value of BF(b — ► S7). As tan/3 increases to 30, 
the HB/FP region stands out as the main green region with low x 2 /dof. Increasing tan/3 
to 45 as in frame a), the HB/FP region remains the most viable, while a yellow corridor of 
neutralino annihilation via s channel A and H appears as splitting the plot. In addition, a 
tiny gray region at low mo and m\/2 has emerged, where m 2 A < 0. As tan/3 increases to 52 
as in frame d.), the m\ < constraint has begun to overwhelm the plot at low mo, pushing 
the Higgs annihilation region to larger parameter values, where a few green points emerge. 
The HB/FP region remains robust. The parameter space becomes completely excluded at 
tan (5 values of 55 and higher. 

The lesson from the < plots is that the most robust region of mSUGRA model 
parameter space is the HB/FP region, where the neutralino has a significant higgsino 
component, so that efficient annihilation (and co-annihilation) of neutralinos can occur in 
the early universe, in spite of quite heavy, multi-TeV values of scalar masses. In fact, these 
scalar masses are sufficiently heavy to suppress possible SUSY CP and flavor violating 
processes, while maintaining naturalness [14]. They thus provide at least a partial solution 



to the SUSY flavor and CP problems. 

If in fact the relic cold dark matter is made of HB/FP neutralinos, can these DM 
particles be detected by direct search experiments? We show the reach of several direct 
detection experiments (CDMS, CDMS2p|/CRESSTp| and Genius |1) for SUSY CDM 



by the black contours, via the spin-independent neutralino-proton scattering rates as cal- 
culated in Ref. |28|]. 3 Similar results are given in Feng, Matchev and Wilczek[fL5|l, although 
their HB/FP region occurs at lower values of mo than ours. 4 It is gratifying to note that the 
most favored regions of parameter space are also accessible to direct search experiments, 



especially large scale experiments such as Genius and Zeplin 4[ 30 1 . 

In Fig. |], we show the same mSUGRA model plane plots, except for \x > 0. In this 
case, the values of BF(b — > 57) and especially a^ USY are more easily accommodated by the 



data] IS]. The frame a.) for tan/3 = 10 is qualitatively similar to the fi < case, with the 
most favorable region again being the HB/FP region. As tan/3 increases to 30, the HB/FP 
region becomes even more promising. In addition, there are some tiny regions along the 
stau co-annihilation border where a low x 2 /dof can be found. The HB/FP region remains 
most promising for tan/3 values of 45 and 52. As in the \i < case, it should be possible 



3 These contours emerge from digitizing the a vs. reach contours presented by the various experi- 
mental groups. We have not scaled the reach contours according to the value of the neutralino relic density. 

4 Our HB/FP region as determined by IS A JET is in good agreement with similar calculations by the 
programs Suspect, SoftSUSY and Sphenop9[. 
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Figure 1: Plot of x 2 /dof for the mSUGRA model in the mo vs. m 1 / 2 plane for /i < 0, Ao = and 
tan/3 = 10, 30, 45 and 52. 



to directly search for the HB/FP dark matter with CDMS2 and especially with Genius 
and/or Zeplin4. 

Finally, we show in Fig. || fi > planes for very large values of tan/3 = 54, 56, 58 and 
60. 5 In frame a.) for tan/3 = 54, the HB/FP region is even more pronounced than in Fig. |2[ 
Meanwhile, a corridor of s-channel Higgs annihilation is opening up at lower values of mo, 
as shown by the green and yellow region. For tan /3 = 56 in frame b.), the m\ < constraint 



5 Such large values of tan/3 fulfill naturalness conditions if one loop corrections are included in evaluating 
the scalar potential fell. 
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Figure 2: Plot of x 2 /dof for the mSUGRA model in the mo vs. m 1 / 2 plane for \i > 0, Ao — and 
tan/3 = 10, 30, 45 and 52. 



has begun to usurp the HB/FP region. In this case, now a broad region of rapid neutralino 
annihilation has opened up where 2m ^ ~ fnA,H- The region intermediate between these 
which is shaded yellow has fi^/i 2 J us t beyond the WMAP constraint. The low mo and 
mi/2 region has a somewhat higher x 2 value due to the value of BF(b — > S7) dropping 
below 2 x 10 -4 . As tan/3 increases to 58 in frame c), a significant region of resonance 
annihilation is evident. Much of it is accessible to direct dark matter search experiments. 
Finally, in frame d.), only a fraction of parameter space remains viable, but none of it with 
a low x 2 /dof value. The entire tjiq vs. m\ii plane is excluded for tan/3 > 62. 
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Figure 3: Plot of x 2 /dof for the mSUGRA model in the mo vs. m 1 / 2 plane for n > 0, Ao — and 
tan/? = 54, 56, 58 and 60. 



To summarize, we have combined the constraints from WMAP on the neutralino relic 
density with constraints from BF(b — ► sj) and (g — 2)^ in a % 2 analysis which determines 
favorable and unfavorable regions of mSUGRA model parameter space. We find the bulk 
neutralino annihilation region (A.) at low mo and m\/2 essentially ruled out by constraints 
from LEP2, BF(b — > sj) and {g — 2)„. In addition, the stau co-annihilation region (B.) 
has only tiny favorable regions, which would require fine-tuning of parameters to satisfy all 
constraints. The HB/FP region (D.) emerges as a significant region satisfying all constraints 
over a wide range of tan (3 values, and also offers at least a partial solution to the SUSY 
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flavor and CP problems. In addition, the neutralino resonance annihilation regions (C.) 
for fj, > and large tan/3 can satisfy all constraints. The favorable HB/FP regions are all 
accessible to direct dark matter search experiments, and much of it should be accessible to 
TeV scale linear colliders [32], since \fi\ is small, and the lightest chargino frequently lighter 
than ~ 500 GeV. The HB/FP region should also be accessible at LHC searches as long as 
7774/2 is not too large, so that gluino pair production occurs at a high enough ratc[33]. 
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